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DYNAMIC PRESSURE LIMITS FOR FLAT PLATES AS RELATED 


TO NUCLEAR FUEL ELEMENTS 
by Roger L. Smith 
Lewis Research Center 

SUMMARY 

Tests with air or helium flowing over single flat plates were performed at room 
temperature to determine the values of dynamic pressure required to cause plate failure. 
Failure was defined as a deflection greater than 0. 010 inch (0. 254 mm). 

The ranges of the dimensions used were as follows: plate thickness, 0. 020 to 0. 040 
inch (0. 508 to 1. 02 mm); plate length, 0. 5 to 4. 0 inches (1. 27 to 10. 2 cm); plate width, 

1. 2 to 2. 0 inches (3. 05 to 5. 08 cm); and flow-channel height, 0. 038 to 0. 091 inch (0. 96 to 
2.31mm). The plate materials selected were lead, aluminum, copper, and steel, which 
provide ranges of Young's modulus of elasticity from 2x10® to 30x10® psi (1. 38x10® to 
20. 7x10® N/cm\ and of Poisson's ratio from 0. 288 to 0. 42. Gas density varied from 
0. 025 to 1. 0 pound per cubic foot (0. 4 to 16 kg/m®). 

The experimental results were compared with an existing theoretical analysis which 
predicts dynamic pressure for sudden, complete plate collapse. The experimental re- 
sults agree with the theoretical prediction for the effect of the material properties, but 
disagree with the prediction for the effect of the plate and flow -channel dimensions. The 
large amount of experimental data scatter and the disagreement of these data with the 
theoretical prediction maybe caused by leading-edge effect, flutter, turbulence, or other 
unknown factors. 


INTRODUCTION 

Various investigators have studied the problem of collapse of parallel -plate fuel ele- 
ments at high coolant velocities in nuclear reactors. Some early experiments on this 
problem are described in references 1 and 2. One explanation for plate collapse was 
proposed by D. R. Miller in reference 3. This analysis postulates that at some critical 
fluid velocity the pressure forces on either side of the plate become unbalanced to the ex- 
tent that collapse occurs. 


In further analytical work, E. B. Johansson (ref. 4) considered the effects of flow 
redistribution and unequal friction pressure drop in the flow channels. Also, J. J. Kane 
(ref. 5) has done some analytical work on the effect of inlet spacing deviations. To con- 
firm Miller’s analysis, W. L. Zabriskie (refs. 6 and 7) conducted experiments with water 
flowing over aluminum and clear plastic plates. These limited experimental results 
agreed reasonably well with the prediction of critical velocities from Miller's analysis. 

This investigation comprises a further and more extensive series of tests. Several 
plate materials and plate dimensions were used. For each combination of material and 
dimensions, several plates were fabricated, and each was individually subjected to in- 
creasing gas flow at room temperature until it failed. Plate failure was defined as either 
a center upstream -edge deflection of 0. 010 inch (0. 254 mm) or complete plate collapse. 
The velocity at failure in terms of dynamic pressure was correlated as a function of ma- 
terial property and of plate geometry parameters suggested by Miller (ref. 3). 

The scope of the experiments included plate thicknesses from 0. 020 to 0. 040 inch 
(0. 508 to 1. 02 mm); plate lengths from 0. 5 to 4. 0 inches (1. 27 to 10. 2 cm); plate widths 
from 1. 2 to 2. 0 inches (3. 05 to 5. 08 cm); flow-channel heights from 0. 038 to 0. 091 inch 
(0. 96 to 2. 31 mm); and gas (helium or air) densities from 0. 025 to 1. 0 pound per cubic 

Q 

foot (0.4 to 16 kg/m ). The plate materials selected were lead, aluminum, copper, and 
steel, which provided ranges of Young's modulus of elasticity from 2x10° to 30x10 psi 
(1. 38x10® to 20. 7x10® N/cm^) and of Poisson's ratio from 0. 288 to 0. 42. 

Edge -effect calculations are discussed in an appendix by A. F. Lietzke. 


SYMBOLS 

A total cross-sectional flow area at the plate 

a plate thickness 

b plate span 

E Young's modulus of elasticity 

g gravitational constant 

h flow-channel height 

l effective plate length, sb 

p pressure 

q dynamic pressure, pV /2g 

2 

q critical dynamic pressure, pV /2g 

V velocity 


2 


V c critical velocity 

w weight flow rate 

p weight density of gas at plate 

v Poisson's ratio 


ANALYTICAL CONSIDERATIONS 

This section describes the basic features and assumptions of Miller's analysis 
(ref. 3) which led to his equation for predicting critical velocity for failure of multiple 
parallel plates and single flat plates. A modification to Miller's equation for critical 
velocity is proposed. This equation is expressed in terms of critical dynamic pressure, 
the parameter of interest herein. 


Miller Analysis 

This analysis "applies the Bernoulli theorem for incompressible flow to determine 
the pressure differences which would be developed across the plates by modification of 
local fluid velocities as a result of plate deflections. Equating these pressure differ- 
ences to those which would be required to produce the corresponding deflections (as in 
a beam) gives an expression for the critical velocity in terms of plate and fluid param- 
eters" (ref. 3). The initial deflections could be caused by gas flow turbulence or im- 
perfections in the hardware; there is always some perturbation present to initiate a 
plate deflection. 

The assumptions were as follows: 

"1. The plates are homogeneous, isotropic, elastic, initially flat . . . , uniform in 
spacing and dimensions, and free of unidentified sources of deformation. 

2. The coolant is incompressible; all channels have the same mass flow; at any 
cross section normal to the longitudinal axis the flow within any channel is uniform; and 
leakage between channels is suppressed. 

3. Plates are broad enough in comparison with their thickness that shear deforma- 
tion is negligible, and are long enough in comparison with their breadth that plates can 
deflect locally without significant redistribution of flow among the coolant channels. 

4. Side plates or supports are rigid. " 

The final equation for the case of multiple parallel flat plates with built-in edges is 
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The analysis indicates that V (and therefore critical dynamic pressure q ) for 

c c 

single flat plates is twice that given by equation (1) for multiple parallel flat plates. The 
equation for a single plate can be written in terms of critical dynamic pressure as 
follows: 
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^The first term on the right side is the plate material parameter (function of material 
■properties) ;Mhe second term is the plate geometry parameter (function of plate and pas- 
sage dimensions). 


Limitations of Miller Analysis 

Nonideal co nditions. - It seems obvious that the assumed ideal conditions cannot 
exist in an actual experiment. The variations from ideal conditions, especially regard- 
ing uniformity of dimensions, materials, and flow, would be expected to cause experi- 
mental plate failure to occur at a lower dynamic pressure than the ideal value predicted 
by the analysis. Since an incompressible fluid was assumed in the analysis, the Mach 
number for the tests reported herein was kept below 0. 4 to minimize compressibility 
effects. 

By assuming the plate to be a wide beam, the analysis neglects any local effects 
occurring in the vicinity of the upstream or downstream edges of the plate. In fact, 
reference 3 states that "If the plates are deflected over a longitudinal distance several 
times the width of the plates, then near the center of the deflected portion of a plate the 
longitudinal variations in deflection can be neglected, . . . .” The analysis then goes on 
to consider the central portion of the plate. 

At the upstream and downstream edges, the plate is no longer, in effect, a wide 
beam, and the deflection at these edges will be greater than that calculated in the anal- 
ysis. The calculations for this greater deflection and its effect on dynamic pressure 
are shown in the appendix. In order to incorporate these effects it is necessary to 
modify equation (2) as given by the following equation from the appendix: 
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Failure criterion . - Since these experiments were intended to measure the dynamic 
pressure required to produce plate failure, some criterion of plate failure was needed. 
Although the analysis indicates that once failure is initiated a sudden and complete 
collapse of the plate would occur, this did not happen in many of the cases reported here 
and elsewhere (ref. 6) . Often the deflection of the upstream edge of the plate would 
gradually increase with an increase in flow. Little deflection occurred in the downstream 
edge. 

A somewhat arbitrary center -upstream -edge deflection of 0.010 inch (0.254 mm) was 
used as the failure criterion for gradual failures to obtain meaningful data for nuclear 
reactor fuel-element design. For this deflection some of the outer fibers of the alumi- 
num, copper, and lead plates will have stretched past their yield point, while no yield 
will have occurred in the steel plates. For the extreme cases less than half the plate 
remains in the elastic condition when deflected 0.010 inch (0.254 mm) at the center. 

Electrical contacts described in Electrical instrumentation section and placed at the 
upstream edge of the plate signaled when the center-upstream -edge deflection had 
reached 0.010 inch (0.254 mm). As mentioned in the INTRODUCTION, earlier tests 
(e. g. , ref. 6) showed that failure nearly always occurred first at this location. 

In the investigation of the effect of upstream -edge support, the plate always failed 
suddenly in a gross manner. This failure was easily observed with the instrumentation, 
so no arbitrary definition of failure was required. 


APPARATUS 

The apparatus used for these tests consisted of the piping system, the plate holder, 
and the instrumentation. 


Piping System 


The piping system is shown in figure 1. The inlet end of the system was connected 
to the laboratory gas supply with 4-inch (10.2-cm) pipe. The outlet end was connected 
to atmosphere outside the building with 6-inch (15.2-cm) pipe. Within this piping sys- 
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tem, the gas flowed through the orifice section, the inlet control valve, the test section 
which contained the plate-holder assembly and test plate, and the outlet control valve. 

Immediately downstream of the inlet connection was a standard ASME orifice section 
for measuring the gas flow rate. Pneumatically operated valves, to control both pres- 
sure level and flow rate, were placed both upstream and downstream of the test plate 
(fig. 1). These 2-inch (5.08-cm) plug-type valves were remotely operated from the con- 
trol panel . 



The plate-holder mounting section was stationed between two pipe tees (fig. 1), each 
of which contained a sight glass to permit plate illumination and observation. This mount-, 
ing section was a spool piece about 1 foot (30. 5 cm) long, with the plate-holder assembly 
mounted to it. For most of the tests, this spool piece was made of steel. However, for 
the tests utilizing high-speed moving pictures of the plate, this spool piece was made of 
clear plastic. 
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Plate-Holder Assembly 


Figure 2 shows the plate -holder assembly. This holder was made of stainless steel 
machined to a thickness of 0. 30 inch (7.62 mm) with bolt holes spaced 1. 75 inches 
(4.44 cm) apart except at the plate location. The plate holder extended 7.75 inches 
(19.7 cm) upstream of the plate for a ratio of length to hydraulic diameter greater than 20. 
This extension allowed ample length to develop undisturbed flow at the upstream edge of 
the plate being tested. The flow-passage surfaces were painted flat black to reduce re- 
flections. 

Two important features of the holder are shown in figure 2(a). One is the contact 
point located 0. 0625 inch (1. 59 mm) downstream from the leading edge of the plate, 
which was used to indicate electrically a plate failure. This contact was also used to 
point support the upstream edge when desired. The other feature is the wall static 
pressure tap, a 0. 0312-inch- (0. 794-mm-) diameter hole, perpendicular to the surface 
and 0. 5 inch (1. 27 cm) upstream of the leading edge of the plate. 



1 -Micrometer screw 



C-67-3177 


(b) Outside view of assembly with micrometer for contact adjustment. 
Figure 2. - Plate holder for single-plate tests. 
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Figure 2(b) shows the outside of the assembly with the micrometer for electrical 
contact adjustment. The other side had a similar adjustable electrical contact, and each 
micrometer head had a setscrew to hold it in position. 

Figure 3 is an enlarged schematic diagram of the plate holder with a typical spacer 
arrangement. The 0.003-inch (0.076-mm) clearance was used so that the plate being 
tested would not be firmly held. Miller's analysis does not allow for the tension forces 
which fixed ends would cause. Variation of the width as well as the thickness of the other 
spacers above and below the test plate produced the desired variations in plate span and 
flow -channel height. Spacers of clear plastic were used when the high-speed moving pic- 
tures were taken. 



Instrumentation 

Four types of instrumentation provided the data required from these tests: pressure, 
temperature, optical, and electrical. 

Pressure instrumentation. - Static pressure taps were placed at four locations 
(fig. 1): two in the orifice flanges, one in the plate-holder assembly, and one in the 
spool piece surrounding the plate holder. These taps were connected by copper tubing to 
indicating devices. Calibrated pressure gages indicated upstream orifice pressure and 
upstream test-plate pressure. Mercury manometers indicated pressure drop across the 
orifice plate and test plate. The pressure gage and manometer indications were recorded 
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manually. The two pressure differentials and the upstream test -plate pressure were also 
measured with pressure transducers and recorded on strip-chart recorder^. 

Tempe rature instrumentation. - Iron-constantan thermocouples were used at two lo- 
cations (fig. 1) to measure orifice air temperature and test-section inlet temperature. 
These temperatures were indicated on a null-type potentiometer and recorded manually. 

Opt ical instrumentation . - The optical system (fig. 1) provided visual monitoring of 
the plate during testing. Light shining through special high-strength glass illuminated 
the test plate. Another glass transmitted the image of plate and contact points to the 
camera by means of the lens and half-silvered mirror. This mirror transmitted part of 
the light to a movie camera and reflected the other part to a television camera. The test 
operator monitored the image on a panel-mounted television screen. 

For the high-speed moving pictures of the side edge of the test plate, the camera was 
placed to view perpendicular to the flow through the clear plastic spacers and spool piece. 
Illumination came from behind so that a silhouette of the plate edge was observed. 

Electrical instrumentation. - Micrometer -adjusted electrical contacts were provided 
in the plate-holder assembly to indicate plate failure. These contacts were connected to 
one of the following indicating devices: an ohmmeter, an oscilloscope, or a buzzer. The 
ohmmeter was used for most of the tests. 

PROCEDURE 

This section describes the procedures used to prepare the plates for testing, to test 
the plates to failure, and to reduce the experimental data to values of dynamic pressure. 


Plate Preparation 

The processing of the plates to be tested was carefully controlled to ensure maximum 
uniformity. The lead, aluminum, and copper plates were rolled to within 0.001 inch 
(0. 0254 mm) of the specified thickness and sheared to size. The aluminum and copper 
plates were then annealed in a furnace while clamped between flat steel plates to make 
them as flat and as uniform in their properties as possible. Each lead plate was flattened 
by placing it on a flat surface and rubbing it. Since annealing takes place in lead at room 
temperature, the lead plate became annealed shortly after flattening. Each steel plate 
was cut from the same piece of flat-sheet spring stock which was the right thickness when 
received; the steel plates were not annealed. Finally, the upstream edge of each plate 
was always deburred but left square. Micrometer calipers were used to measure the 
thickness of each plate. 
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Testing 


The test plate was installed in the plate-holder assembly and the bolts were always 
tightened to a torque of 25 inch- pounds (282 N-cm). The contact points were each succes- 
sively advanced until electrical contact with the plate was made. Each point was then re- 
tracted 0.010 inch (0.254 mm) with the micrometers and was clamped in position with 
setscrews. The holder assembly was then inserted into the test-section spool piece 
(fig. 1). 

The test was started by opening the inlet control valve and setting the desired pres- 
sure level with a small gas flow rate. At this condition a data point was taken for refer- 
ence, and the movie camera was turned on. While constant pressure was maintained, gas 
flow rate was slowly increased until either (1) contact was indicated between the test plate 
and the micrometer -adjusted contact or (2) the plate failed suddenly as indicated visually 
or by system pressure changes. One of these conditions represented plate failure, and 
data were again recorded. The elapsed time from first data point to failure point was 
usually about 2 minutes. 


Data Reduction 

The gas flow rate w was measured by using a standard ASME orifice run. Accord- 
ingly, orifice calibration curves provided the gas flow rate as a function of orifice flow 
coefficients, orifice size, orifice upstream gas density, and measured pressure drop 
across the orifice plate. The temperature-measuring thermocouples were located in 
large-area, low-velocity sections of the pipe; therefore, the temperature measured is 
static temperature. 

With the flow rate known, the velocity was calculated from the continuity equation 


P A 

and the dynamic pressure was calculated from the equation 



The gas density and flow area were evaluated at the test plate. 


( 3 ) 


( 4 ) 
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RESULTS AND DISCUSSION 

The geometry and material effects are presented herein, followed by a comparison of 

the experimental data with the theoretical prediction. Then other effects - plate length, 

gas density, initial plate distortion, and upstream-edge support - are presented. 

Finally, the different modes of plate failure are discussed. 

The data for dynamic pressure at failure are plotted as functions of the variables and 

2 

parameters of equations (2) and (12). The parameter E/(l - v ) is called the plate ma- 

3/4 

terial parameter and is a property of the plate material. The parameter a h/b is 
called the plate geometry parameter since it is a function of the geometry factors . The 
averages presented are arithmetical averages. 


Geometry Effects 

The dimensions which were varied to study geometry effects were plate thickness, 
plate span, and flow-channel height. Aluminum was chosen as the plate material to study 
these variables. The material parameter E/(l - v ) for aluminum has a constant value of 
1. 12x10^ psi (7.73x10^ N/cm^). For these tests, plate length was 2.0 inches (5.08 cm) 
and gas density was about 1. 0 pound per cubic foot (16 kg/m ). 

Plate thickness . - The effect of plate thickness on dynamic pressure at failure is 
shown in figure 4(a). Plate width and flow-channel height were held constant, while plate 
thickness was varied from 0.020 to 0.040 inch (0. 508 to 1.02 mm). The individual data 
points as well as average values for each thickness group are shown. The data show con- 
siderable scatter. The dynamic pressure at failure increased approximately as the sec- 
ond power of the thickness as indicated by the solid line drawn near the mean values. 

Pla te span . - The effect of plate span, or width, on dynamic pressure at failure is 
shown in figure 4(b). The plate span was varied from 1.2 to 2.0 inches (3. 05 to 5. 08 cm) 
with plate thickness and flow-channel height held constant. Again, the experimental data 
show considerable scatter. Dynamic pressure at failure increased approximately as 
l/(b ' ), where b is plate span. 

F low-channel height . - The effect of flow-channel height on dynamic pressure at fail- 
ure is shown in figure 4(c). The heights used were 0.038, 0.064, and 0.091 inch (0.96, 
1.62, and 2.31 mm, respectively), with plate thickness and span held constant. Here 
also, there is considerable scatter in the data. The dynamic pressure at failure in- 
creased approximately with flow-channel height to the 0.61 power. 

Geometry parameter . - The thickness, span, and flow-channel height data, along with 
applicable data points from references 6 and 7, are presented in figure 5. Here the ratio 
of dynamic pressure at failure to plate material parameter is plotted against the plate 


11 


Ratio of dynamic pressure at failure to plate material parameter, q / [ e/( 1 - v 2 )], dimensionless 





r 


60xl0~ 7 Data 



Plate geometry parameter, a^h/b , dimensionless 

Figure 5. - Combination of aluminum data for geometry 
parameter effect. Approximate gas density, 1.0 pound 
per cubic foot (16 kg/m 3 ). 


geometry parameter (see eq. (2)). A solid line is drawn near the average points and has 
a slope of about 0.72 as compared with the theoretical slope of 1.0. 


Material Effects 


The experimental data for four materials, lead, aluminum, copper, and steel, are 

shown in figure 6. Dynamic pressure at failure is plotted against the plate material pa- 

-8 

rameter while the geometry parameter is held constant at 3.21x10" . The solid line 
drawn through the data points shows that the dynamic pressure at failure increases almost 
in direct proportion to the material parameter, as predicted by equation (2). 


Comparison of Experimental Data with Theoretical Prediction 


Figure 7(a) is a combined plot of all the data presented previously for both material 
and geometry effects. The square data points represent the arithmetical mean values for 
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Figure 6. - Combination of data for materials 
parameter effect. Plate geometry parameter, 

3.21X10' 8 ; approximate gas density, 1.0 
pound per cubic foot (16 kg/m 8 ). 

the various materials, while the vertical lines through these points indicate the range of 
data scatter. Also shown is the line of Miller’s theoretical prediction (eq. (2)) and the 
water -flow data points from references 6 and 7, indicated by the circular symbols. The 
solid line drawn through the experimental averages has a slope of about 0.68 compared 
with a slope of 1.0 for the theoretical line. 

The geometry parameter effect is the principal cause of the difference between the 
experimental and theoretically predicted results. A comparison of figures 5 and 6 shows 
that the theoretical analysis predicts the material parameter effect quite well but poorly 
predicts the geometry parameter effect. Figure 7(b) shows the same data as figure 7(a) 
but uses the edge -effect modification applied to the plate geometry parameter. This 
modification is partly a function of plate length l and Poisson's ratio v, neither of 
which were originally included in the geometry parameter. The effect of the correction 
factor resulting from consideration of these edge (and also, therefore, length) effects can 
be seen by comparing the two parts of figure 7. The slope of the experimental line in 
figure 7(b) is about 0. 78 compared with 0. 68 for the line (fig. 7(a)) which was not cor- 
rected for edge effect. Thus, the correction factor helps to give a better theoretical pre- 
diction. Since this prediction still differs from the experimental results, there must re- 
main some unknown factors affecting plate failure. Some of these factors may be other 
leading-edge effects, flutter, and turbulence. These same factors also could contribute 
to the data scatter. 
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Ratio of dynamic pressure at failure to plate material 
parameter, q / [ e/( l - v 2 )], dimensionless 
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(a) Data not corrected for edge effect. (M Data corrected for edge effect. 

Figure 7. - Comparison of theoretical prediction with combined data for all materials and geometries tested. Approximate gas density, 1.0 pound per cubic 
foot (16 kg/m 3 ). 
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Other Effects 


Plate length . - The effect of plate length on dynamic pressure at failure was deter- 
mined by testing plates 0. 5, 2, and 4 inches (1.27, 5.08, and 10.2 cm, respectively) long. 
Also, a 48-inch (122 -cm) length was simulated by placing a flow restrictor at the down- 
stream end of a 4-inch- (10.2-cm-) long plate; the restrictor pressure drop was equiva- 
lent to that for a 48-inch (122 -cm) passage length. Figure 8 shows that measured dynamic 
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0 2 4 6 8 10 122 

Plate length, cm 


Figure 8. - Effect of plate length in flow direction on 
dynamic pressure at failure. Plate material, 
aluminum; plate thickness, 0.020 inch 
(0.508 mm); plate span, 2.0 inches (5.08 cm); 
flow-channel height, 0.064 inch (1.62 mm); 
approximate gas density, 1.0 pound per cubic foot 
(16 kg/m 3 ). The 48-inch (122-cm) length was 
calculated. 
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pressure at failure (divided by a constant material parameter of 1. 12x10 psi (7.73x10 
2 

N/cm )) decreases as plate length increases. These data follow quite closely the theo- 
retical prediction (see appendix) except for short lengths. Flow redistribution may ac- 
count for the divergence of the data at short lengths, since the theory assumes no flow re- 
distribution. 

Gas density . - The density of the hydrogen coolant in a nuclear rocket reactor would 
be orders of magnitude less than the convenient density for room -temperature tests of 

O 

about 1.0 pound per cubic foot (16 kg/m ). Therefore, density-effect tests were under- 
taken. The effect of gas density on dynamic pressure at failure was evaluated using heli- 
um and air at several pressure levels. The result is shown in figure 9. The influence of 
gas density appears relatively minor compared with other unknown factors that cause 
great data scatter. 
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Figure 9. - Effect of gas density on dynamic pressure at failure. Plate 
material, aluminum; plate thickness, 0.020 inch (0.508 mm); plate 
span and length in flow direction, 2.0 inches (5.08 cm); flow-channel 
height, 0.094 inch (1.62 mm). 
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Figure 10. - Effect of supported upstream edge. Plate length, 
4.0 inches (10.2 cm); approximate gas density, 1.0 pound 
per cubic foot (16 kg/m 3 ). 



Initial plate dist ortion. - The initial micrometer settings of the contact points were 
tabulated for the data from five different geometry parameters . From this data the off- 
center position of the plate was calculated. In each group the most off-center plate did 
not give a data point at the extremities of the scatter. Thus, initial lack of flatness is 
apparently not an important factor in the scatter of the data points. 

Upstream -edge supp ort. - In an attempt to eliminate leading-edge effect, seven se- 
ries of tests were made with the upstream edge supported at the center by the electrical 
contacts. These contacts were about 0.03 inch (0.762 mm) in diameter, and the plates 
were 4 inches (10.2 cm) long. Each series consisted of either five or six individual tests. 

Figure 10 shows this data, along with the theoretically predicted line, where the span 
was not considered to be changed by the support. The slope of a line through this data is 
about the same as that for a line through the data taken with no upstream -edge support 
(fig. 7, p. 15). Thus, the effects of materials and geometries are about the same re- 
gardless of upstream-edge support. Upstream-edge support increased the dynamic pres- 
sure at failure by a factor of about 2. 4. Scatter was considerably reduced, and all fail- 
ures were sudden as indicated by visual or pressure indications. Thus, the condition of 
the leading edge is important, and therefore the theoretical prediction (eq. (2)), which 
neglects the edge condition, does not accurately predict the dynamic pressure at failure. 


Modes of Plate Failure 

Two basic modes of plate failure were observed during tests without upstream-edge 
support. The first was a sudden and complete failure of the test plate from its neutral 
position. The second mode of plate failure was a gradual deflection of the plate from its 
neutral position with increasing flow, until it touched the electrical contact or, on occa- 
sion, developed severe vibration, which produced intermittent electrical contact followed 
by complete failure. 

For the aluminum plates tested without upstream support, about half the failures 
were sudden. For the copper and lead plates tested, also without upstream support, most 
failures (88 percent) were gradual. Although the data scatter was large for either type 
failure, the average dynamic pressure for gradual failures of aluminum plates was only 
about 5 percent lower than that for sudden failures. Since the difference is so small, 
both sudden and gradual failures are discussed together. 

During the tests with upstream-edge point support, the failures were always sudden. 
In an attempt to determine the mode of sudden failure initiation, high-speed moving pic- 
tures were taken during three of these tests. Figure 11 shows a sequence of high-speed 
film frames taken during initiation of failure of an aluminum plate. The pins supporting 
the upstream edge are on the right side of each picture; the arrow shows the flow direc- 
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Figure 11. - High speed motion-picture frames of plate failure. 


tion. Only about 2 inches (5. 08 cm) of the upstream portion of the plate is shown. The 
film frame rate was in the range of 8000 to 9000 frames per second. 

In the figure, the initial deflection appears to be downward at the upstream edge of 
the plate (right side of the figure) and upward near the center of the plate. Following this 
a wavelike motion of the plate is observed, moving downstream along the plate at a veloc- 
ity of about 450 to 500 feet per second (137 to 152 m/sec) as indicated by the film frame 
rate. This velocity range is the same as that of the air passing over the plate. The 
complete sequence of film frames shows less than one cycle of plate motion before the 
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plate comes to rest in the flowing air. The other two similar tests, not illustrated, 
showed as many as five complete cycles before the motion stopped. Thus, what appeared 
to be a sudden failure was really something like flutter followed by permanent deforma- 
tion. 


SUMMARY OF RESULTS 

The results of an experimental investigation to determine the dynamic pressure at 
failure for single flat plates are presented. These results were compared with an exist- 
ing theoretical analysis which predicts critical dynamic pressure (the dynamic pressure 
which should cause sudden, complete plate collapse). The experimental results agree, 
in general, with the theoretical prediction for the effect of the material properties, but 
disagree with the prediction for the effect of the plate and flow-channel dimensions. A 
modification based on edge effect was made to the theoretical analysis which improved 
the accuracy of the prediction. Also, a large amount of data scatter was present in the 
experimental results. Some probable reasons for the lack of agreement with the theo- 
retical prediction and for the data scatter are leading-edge effect, flutter, and turbu- 
lence, as well as other factors presently unknown. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 26, 1967, 

122-28-02-04-22. 
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APPENDIX - EDGE-EFFECT CALCULATIONS 
by A. F. Lietzke 


Since Miller considers a central portion of the plate, his equation does not account 
for what happens at the upstream and downstream edges. There are two effects which 
should be considered. The first is the increased deflection of the plate (as a beam) near 
the edges. The second is the angle of attack produced by the deformation of the edge of 
the plate when the whole plate bends as a beam. 

An estimate of the edge effect is obtained by using an equation from reference 8. 

This equation describes the position of any point along the surface of a prismatical bar in 
pure bending. Thus, the position and angle of attack of the edge of the beam can be calcu- 
lated and their effect on pressure determined. 


Area Change Modification 

This calculation was performed to obtain the ratio by which the flow -area change at 
the upstream edge is altered by the edge effect. To do this requires knowledge of the po- 
sition of the edge relative to the average position of the plate. 

Figure 12 shows a prismatical bar bent by two equal and opposite couples. This 
situation is considered to be analogous to the bending of a beam with built-in ends. Ref- 


Flow 



Figure 12. - Prismatical bar bent by two equal and opposite couples. 


erence 8 gives an equation for the position of any point on the surface of the bar. The 
following three qualifications are made in using the equation: 

(1) The shape of the cross section is independent of the position along the bar in the 
z -direction. 

(2) The equation describes the centerline of the plate as shown in figure 13. 

(3) The l dimension cannot exceed the b dimension, where b is the plate span 
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Figure 13. - Enlarged view of prismatical bar showing relative deflections. 


perpendicular to l, the length in the flow direction. Reference 9 treats the buckling of 
flat sheets with sides restrained and concludes that when the length is greater than the 
span, the plate ’’tends to act as a series of square plates with buckles in each square, 
since it is prevented from buckling in a single long wave. ” Thus, the plate tends to 
assume a wavelike shape along its lengthwise centerline with nodes separated by the dis- 
tance b; therefore, l can never be greater than b in the equations about to be devel- 
oped. 

When applied to the situation shown in figure 13, the equation of this parabolic curva- 
ture is 


2R 

where x and y are the coordinates of any point on the surface using the axes shown in 
figure 12 and R is the radius of curvature of the bar after bending. Referring to fig- 
ure 13, the position of the end Xq will be obtained when this equation is evaluated at 
y = l / 2 which yields 


x 


0 ~ 



The average position of the centerline x„_, is determined by 

civ 
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X 


av 


Width = 


Area 

Length 




vl 

24R 


The absolute value of the difference 


. _ x I = Hi 2 - Hi 2 = Hi 2 

av 01 24R 8R 12R 


The ratio of flow -channel heights caused by this deflection is 

D* D + l x 0- x av| _ 1 + 

D D 12RD 

where D is the average deflection of the plate centerline and D' is the position of the 

edge of the plate. From beam theory, R = (EI/M) and D = (p^b^/384EI) where I is the 

moment of inertia of the plate cross section (in. ^ (cm^)) and M is the bending moment 

2 2 

(in. -lb (N-cm)). In terms of the moment at midspan, where M = (pfb /24), D = (24Mb / 
384EI) and then 


IT = j | vl 2 M(384EI) _ j + Avl 2 
D 12EI(24Mb 2 ) 3b 2 


Therefore, 


AA’ = 1 + 4ul 2 
AA 3b 2 


( 6 ) 


where AA is the area change corresponding to the average deflection D of the plate and 
A A’ is the area change at the edge of the plate. 


Lift Modification 

To obtain the angle of attack a introduced by this edge effect, start with equa- 
tion (5) 
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( 5 ) 



2R 

Differentiating to get the slope at the edge gives 

<*= JL 2 y 

dy 2R 


where y= (i/2) (fig. 12) so that 

dx _ vl 
dy 2R 

But R = (EI/M) and 

M = N * 2 | Nj 2 Ni 2 _ Nj 2 _ pib 2 
12 4 8 24 24 

where N is the beam load in pounds per foot (N/cm). Since the angle is small, 

a *tan a = *i = n = ^pib 2 = pi 2 b 2 (r 

dy 2R 2EI(24) 48EI 

From the deflection equation, D = (pibV384EI). Rearranging yields El = (pibV384D). 
Substituting this expression into equation (7) results in 

a _ ^pi 2 b 2 (384D) _ 8^jD (8) 

48pi b 4 b 2 

But D = (AA/b), where A = hb; therefore, b = (A/h) and D = (AAh/A). Substituting this 
expression into equation (8) gives 


a = 


Qvl AAh 
b 2 A 


The equation for the lift pressure on a wing for a small angle of attack is 


( 9 ) 
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P =ex! 2 , a 

2g 


Substituting for or from equation (9) gives 

p _pV 2 ff 8Wh AA 
g b 2 A 


( 10 ) 


Combination of Modifications 

The equation for the pressure force caused by a small deflection of the plate is 
(ref. 3) 


p= 2£V^aa (11) 

g A 

where AA is the flow -area change caused by the deflection, and A is the flow area be- 
fore the change. This pressure is increased by the lift force as follows (see eqs. (10) 
and (11)): 


_ 2 pV 2 AA . pV 2 „ 8ul h AA 

p total ~ — 7~ + ~~ 77 T 

g A g v.2 A 


2pV z AA L + 4jr£l h 


Ag 


Including the area-change correction from equation (6) gives 


Ag 


p _ 2 pV z AA L 1 4n vl + 4 vl 


3b 


When this expression is put into the equation for critical dynamic pressure, there results 
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( 12 ) 



remembering that l cannot exceed b. 
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